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Abstract. We have previously presented 5.7 a evidence for a smaller a at redshifts 
0.2 < Zabs < 3.7 from a sample of 128 Keck/HIRES quasar absorption systems: 
Aa/a — (— 0.57±0.10) x 10~^. A non-zero /^a/a manifests itself as a distinct pattern 
of shifts in the measured absorption line wavelengths. The statistical error is now 
small: we do detect small line shifts in the HIRES data. Whether these shifts are due 
to systematic errors or due to real variation in a is now the central question. Here 
we summarize the two potentially most important systematic effects: atmospheric 
dispersion and isotopic abundance evolution. Previously, these have been difficult 
to quantify/model but here we find that neither of them can explain our results. 
Furthermore, the HIRES spectra themselves contain no evidence for these effects. 
Independent measurements of Aa/a with a different telescope and spectrograph are 
now crucial if we are to rule out or confirm the present evidence for a variable a. 
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1. Introduction 



In Murphy et al. (2002) and Webb et al. (2003) we presented new evi- 
dence for a smaller fine structure constant, a = /he, in high redshift 
quasar absorption clouds: /S.a/a = (—0.57 it 0.10) x 10"^. This prelimi- 
nary result is the weighted mean from a many-multiplet (MM) analysis 
of 128 absorption systems and represents a significant improvement on 
our previously published result from 49 systems (Murphy et al., 2001a; 
Webb et al., 2001): Aa/a = (-0.72 ± 0.18) x 10"^. Such a potentially 
fundamental result requires extreme scrutiny, including a thorough con- 
sideration of possible systematic errors. In Murphy et al. (2001b) we 
explored a wide range of possibilities: laboratory wavelength errors, 
wavelength calibration errors, atmospheric dispersion effects, uniden- 
tified interloping transitions, isotopic ratio and/or hyperfine structure 
effects, intrinsic instrumental profile variations, spectrograph tempera- 
ture variations, heliocentric velocity corrections, kinematic effects and 
large scale magnetic fields. All but two of these - atmospheric dispersion 
effects and isotopic abundance variation - were shown to be negligible. 
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Here we present an updated analysis of these two effects. We also 
present new results from a line removal test which constrains systematic 
errors associated with individual transitions or species. 

2. Atmospheric dispersion effects 

Fig. |l|a illustrates how atmospheric dispersion results in an angular sep- 
aration between different wavelengths as they enter the spectrograph 
slit. This separation increases with increasing zenith angle, ^. For the 
particular optical design of the Keck/HIRES, the component of angular 
separation across the slit (i.e. in the wavelength direction as projected 
from the CCD) will cause a compression of the quasar spectrum. Such 
a compression will mimic Aa/a < for the Mg/Fe il systems observed 
at Zabs < 1-8, the effect being complicated at high-z by the diversity of 
q coefficients^ (see Fig. ^). A wavelength dependent asymmetry will 
also be introduced into the instrumental profile (IP) due to truncation 
of the seeing discs by the slit edges. These effects will not be applied 
to the thorium-argon (ThAr) calibration spectra since the ThAr lamp 
illuminates the slit uniformly in all cases. 

An image rotator is usually employed to hold the slit axis per- 
pendicular to the horizon as the quasar is tracked across the sky. 
This eliminates the compression and IP asymmetry effects. Only 51 
of our 128 absorption systems were observed in this way; 77 systems 
were observed before August 1996 when Keck/HIRES was fitted with 
an image rotator. However, the values of Aa/a for the 77 affected 
and 51 unaffected systems are consistent, with weighted mean values 
Aa/a = (-0.54±0.14)xl0-5 and (-0.62±0.15) x 10"^ respectively. We 
also expect to observe a correlation between Aa/a and ^, especially for 
the low-z Mg/Fe ll systems (see Fig. Qa). This expectation is confirmed 
in the simulations below. However, no such correlation is observed in 
our quasar data. Therefore, the data themselves suggest a negligible 
effect on Aa/a due to atmospheric dispersion effects. 

Nevertheless, we used the atmospheric and telescope parameters for 
each quasar observation to model the potential effect of atmospheric 
dispersion on Aa/a. The angular separation of two wavelengths enter- 
ing the slit is determined by the refractive index of air at the telescope 
(which depends on temperature, pressure and relative humidity). The 
slit angle with respect to the vertical, 9, is equal to for the pre-rotator 
Keck/HIRES. We applied the resulting wavelength shifts to a high S/N 
simulation of our data and determined Aa/a, giving a correction which 

^ The q coefficient measures the sensitivity of a transition wavelength to changes 
in Q. See (Dzuba et al., 2002) for details and updated values. 
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Figure 1. Summary of known, potentially significant systematic errors, (a) Earth's 
atmosphere disperses seeing discs of different wavelength along the wavelength direc- 
tion (A) of the spectrograph slit, potentially compressing the quasar spectrum. The 
seeing discs are also truncated by the slit edges, leading to wavelength dependent 
IP asymmetry, (b) Simulated Mg I A2852 absorption line (histogram) with only the 
^*Mg isotope. A fit to the data (solid line) with the full isotopic structure (tick-marks; 
separations exaggerated by a factor of 25), using the terrestrial isotopic ratios, 
induces an artificial line shift, (c) The upper panel corrects Aa/a for atmospheric 
dispersion effects and the middle panel corrects for strong isotopic ratio evolution 
of ^'^'^^'Mg and ^^'^°Si. Both corrections are included in the lower panel. 



we apply to the values measured from the quasar spectra. We compare 
the corrected values of Aa/a with the uncorrected values in Fig. |l|c. 
Overall, our model of atmospheric dispersion has an opposite effect on 
the low-z systems compared to the high-z systems: it can not explain 
the observed deviation from Aa/a = 0. Once we apply the correction, 
the apparent evolution of a with look-back time is enhanced. 



3. Isotopic ratio evolution 

Mg, Si, Cr, Fe, Ni and Zn have several stable isotopes and so the 
transitions used in the MM method will have unresolved isotopic struc- 
ture. However, only the isotopic structures of the Mg lines have been 
determined experimentally (Hallstadius, 1979; Drullinger et al., 1980). 
The Mg and Si transitions involve similar energy levels and so, in Mur- 
phy et al. (2001b), we estimated their isotopic structure by scaling the 
Mgii A2796 structure by the mass isotopic shift. We fitted the Mg and 
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Si absorption lines with these isotopic structures using the terrestrial 
values of the isotopic abundances (Rosman and Taylor, 1998). However, 
Galactic observations of Mg (Gay and Lambert, 2000) and theoretical 
models of Si (Timmes and Clayton, 1996) in stars strongly suggest that 
only the ^^Mg and ^^Si isotopes will exist in quasar absorption clouds 
with significant abundances. Fig. |l]b illustrates the effect this would 
have on a measured Mgi A2852 absorption line position (the isotopic 
separations are exaggerated by a factor of 25). 

To obtain an upper limit on the effect of isotopic ratio evolution 
on Aa/a, we refit the Mg and Si transitions using only the ^^Mg 
and ^*Si isotopes, i.e. zero abundances of ^^'^^Mg and ^s.sogj ^j^g 
absorption clouds. We compare these new, corrected values of Aa/a 
with the uncorrected values in Fig. ||c (middle panel). As expected, 
the potential effect is most prominent in the low-z Mg/Fe ii systems. 
Note that Aa/a is more negative after the correction: isotopic ratio 
evolution can not explain our results. Although the isotopic structures 
of the other transitions in our analysis are not known, they should be 
more compact since the mass isotopic shift is inversely proportional 
to the square of atomic mass. The line removal test described below 
largely rules out such strong effects. 



4. Line removal tests 

A non-zero Aa/a manifests itself as a distinct pattern of line shifts 
(Fig. ^). Thus, with a large set of fitted transitions, one or more 
transitions can be removed from a given system and Aa/a will remain 
well-constrained. In the absence of systematic errors associated with 
the transition removed or with the wavelength scale of the quasar 
spectra, Aa/a should not change systematically for all systems. For 
example, removing a single transition (e.g. Mgii A2796) from our fits 
to the quasar data directly probes possible systematic errors due to 
line blending, isotopic/hyperfine structure effects and errors in the 
laboratory wavelengths. Removing entire species (e.g. all transitions 
of Sin) also tests for isotopic/hyperfine structure effects. 

Given a set of fitted transitions, we only remove transitions that can 
be removed. To clarify this, consider removing Mgii A2796 (Fig. |2|b). 
For the low-z Mg/Fe ii systems, a well-constrained value of Aa/a can 
only be obtained when at least one Mg line is fitted since the difference 
between the q coefficients for the various Fe ii transitions is small (see 
Fig. ^). If the only available anchor line is Mgii A2796 then we can 
not remove it from the system. Similarly, we only remove the entire 
Mgii species if the Mgi A2852 line is fitted. 
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Figure 2. (a) Distribution of q coefficients for all transitions in the MM method. 
The simple arrangement for the low-z Mg/Fe ll systems indicates susceptibility to 
low-order distortions of the wavelength scale, (b) Removing the Mg I A2796 line is 
possible for n = 61 systems. The effect of removal on these systems (upper panel) and 
on the entire sample (lower panel) is small, indicating negligible systematic errors 
associated with Mgl A2796. (c) Full line removal results. The transition/species 
removed is given on the vertical axis along with the number of systems, n, where 
removal was possible. The left panel compares the weighted mean value of Aa/a 
before (dotted error bar) and after (solid error bar) removal for these n systems. The 
right panel shows the effect of removal on Aa/a for the entire sample, i.e. including 
systems where the transition/species could not be removed. 



In Fig. ^ we present line removal results for all transitions and 
species. The left panel shows the overall effect line removal has on the 
n systems where line removal is possible (cf. upper panel of Fig. |2|b). 
The right panel shows the effect on the entire sample (cf. lower panel 
of Fig. |2|b). We see no evidence for systematic errors associated with 
particular transitions or species. 



5. Conclusions 



We have extended our previous analysis of systematic errors (Murphy 
et al., 2001b) to include the new, significantly augmented sample of 128 
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absorption systems outlined in Murphy et al. (2002) and Webb et al. 
(2003). We cannot explain the observed value of Aa/a = (—0.57 ± 
0.10) X 10~^ with known systematic errors. Here we have outlined our 
analysis of the two most important of these: atmospheric dispersion 
effects and possible isotopic ratio evolution. The former is an instrumen- 
tal effect that should have affected 77 of the 128 absorption systems. 
However, the data themselves suggest that atmospheric dispersion had 
a negligible effect on Aa/a. Furthermore, modelling indicates that 
the resulting low-order distortion of the wavelength scale should have 
shifted the low- and high-z values of Aa/a in opposite senses, clearly 
at odds with our observations. An upper limit on the effect on Aa/a of 
strong isotopic ratio evolution was obtained by assuming zero isotopic 
abundances for ^^'^^Mg and ^^''^''Si in the quasar absorption clouds. 
Although the isotopic structures of the other transitions used in the 
MM method are not known, strong isotopic ratio evolution for these 
elements is strongly constrained by the line removal test (Fig. |2|c). 

All 128 absorption systems were observed with the Keck/HIRES. 
Although we have ruled out many instrumental systematic errors, in- 
dependent measurements of Aa/a must be carried out with a different 
telescope and spectrograph. Using the MM method to analyse VLT and 
Subaru spectra will be a crucial check on possible systematic errors. 
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